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1.  Introduction 


The  CTH  shock  physics  code,  developed  by  Sandia  National  Laboratories,  is  widely 
used  for  defense  and  industrial  applications  involving  shock  waves,  large  material 
deformations,  and  material  failure.1  In  it’s  original  configuration,  CTH  employs  a 
2-step  solution  method;  a  Lagrangian  step  followed  by  a  remap  step.  The  remap 
step,  also  referred  to  as  the  advection  step,  uses  operator-splitting  techniques  to 
replace  multidimensional  equations  with  a  set  of  1 -dimensional  equations.  High- 
resolution  interface  trackers  are  used  to  minimize  material  dispersion  resulting  from 
advection.2 

Inherent  limitations  in  computational  accuracy  resulting  from  advection  have  pro¬ 
vided  the  motivation  for  implementing  a  new  computational  method  into  CTH.3 
As  a  result,  the  Lagrangian  marker  method  has  recently  been  implemented  into 
CTH  and  offers  an  alternative  to  the  legacy  Eulerian  method  discussed  in  the  pre¬ 
vious  paragraph.  The  primary  advantage  of  this  technique  over  the  legacy  Eulerian 
method  is  the  application  of  a  Lagrangian  numerical  approach  to  the  motion  of 
marker  fields  rather  than  using  advection  to  move  material  through  the  mesh.  The 
intended  result  is  a  computational  method  with  the  advantages  of  a  traditional  Eu¬ 
lerian  method — large  material  deformations  and  failure — without  the  deficiencies 
associated  with  material  advection.4 

The  implementation  of  the  Lagrangian  marker  method  into  CTH  marks  a  significant 
change  to  the  functionality  of  the  code.  The  introduction  of  the  new  computational 
method  requires  evaluation  before  it  can  be  adopted  as  a  replacement  for  the  legacy 
Eulerian  solution  method.  The  US  Army  Research  Laboratory  has  begun  an  ini¬ 
tiative  to  evaluate  the  CTH  Lagrangian  marker  method  for  problems  in  terminal 
ballistics.  The  approach  taken  is  to  begin  with  fundamental  problems  in  ballistics, 
with  trusted  experimental  results,  and  verify  the  proper  solution  to  these  problems 
before  advancing  to  the  next  set  of  more  complicated  ballistics  problems.3 

This  report  describes  the  initial  step  in  the  evaluation  of  the  CTH  Lagrangian  marker 
method.  The  Taylor  impact  problem5  was  used  as  the  basis  of  the  study  because  of 
the  availability  of  high-quality  experimental  data  and  the  sensitivity  of  the  experi¬ 
mental  result  to  the  constitutive  behavior  of  the  impactor,  providing  an  opportunity 
to  verify  the  proper  implementation  of  constitutive  models  for  metals  in  the  newly 
implemented  method. 


1 


The  computational  study  involved  the  execution  of  3  different  sets  of  numerical  sim¬ 
ulations  to  match  selected  Taylor  impact  experiments  previously  reported  in  litera¬ 
ture.  CTH  simulations  were  performed  using  the  legacy  Eulerian  method  and  then 
matching  simulations  were  performed  using  the  Lagrangian  marker  method.  These 
pairs  of  simulations  served  to  isolate  the  influence  of  the  marker  method  on  the  solu¬ 
tion  for  each  case  considered.  Finally,  a  set  of  Arbitrary  Lagrangian/Eulerian  (ALE) 
simulations  was  performed  for  comparison  using  the  arbitrary  Lagrangian/  Eulerian 
code,  ALE3D,  developed  by  the  Lawrence  Livermore  National  Laboratory.6 

2.  Taylor  Impact  Experiments 

The  Taylor  impact  experiments  used  in  the  study  involved  metallic  right-circular 
cylinders  fired  into  a  nearly  rigid  surface.  The  impactors  considered  for  the  study 
were  made  of  oxygen-free,  high-conductivity  (OFHC)  copper  and  4340  steel  (Rc30 
hardness)  with  impact  velocities  ranging  from  130  to  343  m/s. 7,8  Comparing  the 
initial  and  final  shapes  of  the  impactors  provides  valuable  information  on  the  high- 
rate  constitutive  behavior  of  the  materials.  Table  1  provides  a  summary  of  the  ex¬ 
perimental  configurations  used  in  the  study  and  the  final  impactor  shapes  from  the 
experiments  selected.  The  values  listed  in  the  table  are  the  striking  velocity  (Vs), 
the  initial  length  (L0)  and  diameter  (D0)  of  the  impactor,  the  final  length  (L)  and  di¬ 
ameter  (D)  of  the  impactor,  and  the  final  “bulge”  of  the  impactor  (W).  The  bulge  is 
defined  by  Johnson  and  Holmquist  as  the  postevent  diameter  of  the  impactor  mea¬ 
sured  from  a  distance  of  20%  of  the  initial  length  (0.2L0)  from  the  impact  face  of 
the  body.9 


Table  1  Summary  of  Taylor  impact  experiments7,8 


Material 

Vs 

(m/s) 

L0 

(cm) 

Do 

(cm) 

L 

(cm) 

D 

(cm) 

W 

(cm) 

130 

2.54 

0.762 

1.956 

0.991 

0.899 

OFHC  Copper 

190 

2.54 

0.762 

1.621 

1.351 

1.007 

208 

2.54 

0.762 

2.253 

1.062 

0.831 

4340  Steel.  Rc30 

343 

0.81 

0.762 

0.631 

1.147 

0.912 

In  a  previous  study,  Johnson  and  Holmquist  developed  a  metric  to  provide  a  single 
quantitative  comparison  between  Taylor  impact  experiments  and  corresponding  nu¬ 
merical  simulations.  This  metric  is  the  average  error  of  the  final  length,  diameter, 
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and  bulge  between  the  experiment  and  simulation  and  is  represented  in  Eq.  I,9  in 
which  AL,  AD,  and  AW  are  the  differences  between  the  experimental  and  com¬ 
putational  results.  In  the  current  study,  this  metric  is  used  as  the  primary  means 
of  evaluating  the  CTH  Lagrangian  marker  method  simulation  results  against  the 
experimental  data  and  the  other  computational  methods. 

1  /  |AL|  \AD\  \AW\ 

3  V  L  D  W 

3.  Constitutive  and  Fracture  Models 

Three  commonly  used  metal  constitutive  models  were  evaluated  in  the  study  to 
verify  the  proper  implementation  of  each  into  the  CTH  Lagrangian  marker  method. 


3.1  Johnson-Cook 

The  Johnson-Cook  constutive  model,  shown  in  Eqs.  2  and  3  provides  a  representa¬ 
tion  of  the  dynamic  yield  strength  of  metals  as  a  function  of  equivalent  plastic  strain 
(ep),  plastic  strain  rate  ( ip ),  and  temperature  (T).7  The  model  parameters  (A,  B,  n, 
C,  Tmpit,  and  m)  for  OFHC  copper  and  4340  steel  (Rc30)  used  in  the  study  were 
derived  by  Johnson  and  Cook.7 

a  =  {A  +  Be;) 

T  = 


1  +  Cln—  )  (1  -  T'm^ 
£o. 


(2) 


T  -  T 

J-  -L 


Tmelt,  Tr00rn 


(3) 


3.2  Zerilli-Armstrong 

The  Zerilli-Armstrong  constitutive  model  is  a  dislocation-mechanics-based  model 
that  incorporates  the  effects  of  work  hardening,  strain-rate  hardening,  and  ther¬ 
mal  softening  based  on  thermal  activation.10  The  model  takes  into  consideration 
grain  orientation  and  includes  forms  for  both  face-centered-cubic  (fee)  and  body- 
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centered-cubic  (bcc)  metals  and  takes  the  general  form  of  Eq.  4,  in  which  ep  is  the 
equivalent  plastic  strain,  ep  is  the  plastic  strain  rate,  and  T  is  the  temperature. 


a  =  A  +  (Cl  +  C2^£p)  e(-C3T+C4T  lnep)  +  c&£n  _  (4) 

The  model  can  be  configured  for  fee  and  bcc  metals  by  the  selection  of  the  model 
parameters.  For  fee  metals,  the  Ci  and  C5  parameters  are  set  to  zero  resulting  in 
the  form  of  Eq.  5  and  for  bcc  metals,  the  C2  parameter  is  set  to  zero  resulting 
in  the  form  of  Eq.  6. 10  The  Zerilli-Armstrong  model  parameters  for  OFHC  copper 
were  obtained  from  Zerilli  and  Armstrong,10  and  the  parameters  for  4340  steel  were 
obtained  from  Gray  et  al., 11  with  the  Ci  parameter  adjusted  to  obtain  the  appropriate 
quasi-static  yield  strength  for  a  hardness  of  Rc30. 


afcc  =  A  +  C2^e(-C3T+C4Tln£p)  .  (5) 

abcc  =  A  +  Cle(~C3T+C4Tln)  +  C5en  .  (6) 


3.3  Steinberg-Lund _ 

The  Steinberg-Lund  constitutive  model12  is  a  strain-rate-dependent  model  derived 
from  the  Steinberg-Guinan  high- strain-rate  constitutive  model.13  The  Steinberg- 
Lund  model  extends  the  validity  of  the  original  Steinberg-Guinan  model  to  strain 
rates  as  low  as  10_4/s.  In  this  model,  the  yield  strength  is  defined  as  a  function 
of  plastic  strain  (ep),  plastic  strain  rate  (e'p),  pressure  (P),  and  temperature  (T)  as 
described  in  Eqs.  7-10. 


a  —  [< 7t{sp ,  T )  +  oa(  1  +  P(£p  +  g))"] 


~G(P,T)- 

G0 


(7) 


T) 


(8) 
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B(T  -  Tr) 


(9) 


G(P,T) 

G0 


1  + 


AP 

T~ 
T]  3 


n  = 


p_ 

Po 


(10) 


The  yield  strength  in  Eq.  7  is  expressed  in  terms  of  the  thermally  activated  com¬ 
ponent  of  the  yield  strength  (crT,  which  is  a  function  of  the  plastic  strain  rate  and 
temperature),  the  thermally  independent  component  of  the  yield  strength  {a a),  and 
the  local  shear  modulus  (G),  which  is  a  function  of  the  pressure  and  temperature. 
The  thermally  activated  yield  strength  is  obtained  from  Eq.  8,  which  is  an  expres¬ 
sion  of  the  plastic  strain  rate  as  a  function  of  the  thermally  activated  yield  strength 
and  temperature.  The  ratio  of  the  local  shear  modulus  to  the  reference  shear  modu¬ 
lus  is  given  in  Eq.  9,  in  which  q  is  the  material  compression  as  expressed  in  Eq.  10. 
The  terms  £,  and  Tr  represent  the  initial  plastic  strain  and  room  temperature,  re¬ 
spectively.  The  remaining  terms  are  model  constants  obtained  from  Steinberg.14 


3.4  Johnson-Cook  Fracture  Model 

The  Johnson-Cook  fracture  model8  was  employed  in  all  the  CTH  Euler  and  Marker 
simulations  in  the  study  and  in  the  ALE3D  simulations  using  the  Johnson-Cook  and 
Zerilli- Armstrong  constitutive  models.  The  model  defines  the  local  failure  strain 
(c  j)  as  a  function  of  the  pressure-stress  ratio  (<r*),  the  normalized  strain  rate  (i*), 
and  the  normalized  temperature  (T')  as  shown  in  Eq.  11,  in  which  the  normalized 
temperature  is  the  same  as  that  defined  in  Eq.  3.  The  values  D,  through  D5  are 
material  model  constants  derived  by  Johnson  and  Cook.7 


£f  =  {p i  +  D2eD3rT*)  (1  +  DAln{e*))  (1  +  DbT')  .  (11) 

Damage  in  the  Johnson-Cook  fracture  model  is  defined  as  the  ratio  of  accumulated 
strain  to  the  failure  strain  as  defined  in  Eq.  1 1 .  When  sufficient  strain  has  accumu¬ 
lated  such  that  the  damage  reaches  a  value  of  1.0,  the  material  at  that  location  is 
considered  to  be  failed  and  can  no  longer  support  shear. 
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4.  Description  of  Numerical  Simulations 


All  simulations  performed  in  the  study  were  2-dimensional  (2-D)  with  axial  sym¬ 
metry.  The  CTH  Euler  simulations  employed  a  fixed  mesh  with  a  uniform  cell  size 
of  0.381  mm,  resulting  in  10  cells  across  the  radius  of  the  impactor.  In  the  CTH  La- 
grangian  marker  method  implementation,  the  markers  use  the  CTH  Eulerian  mesh 
as  a  background  grid  to  track  marker  motion.4  As  a  result,  the  resolution  of  a  marker 
simulation  is  controlled  by  both  the  Eulerian  mesh  definition  and  the  density  of 
markers  per  mesh  cell.  In  the  CTH  marker  implementation,  marker  density  is  de¬ 
fined  as  the  number  of  markers  per  linear  dimension  in  a  cell.  Thus,  for  a  2-D 
simulation,  the  number  of  markers  per  cell  is  the  square  of  the  number  of  markers 
per  linear  dimension,  and  in  a  3 -dimensional  (3-D)  simulation,  the  number  of  mark¬ 
ers  is  equal  to  the  cube  of  the  number  of  markers  per  linear  dimension.  From  this 
logic,  it  is  clear  that  the  computational  cost  of  a  marker  simulation  can  be  strongly 
dependent  upon  the  marker  density. 

For  this  study,  all  CTH  marker  simulations  used  the  same  mesh  as  the  Eulerian  sim¬ 
ulations.  An  initial  set  of  marker  simulations  was  performed  using  3  and  6  markers 
per  linear  cell  dimension,  resulting  in  9  and  36  total  markers  per  cell  in  the  2-D 
simulations.  This  approach  was  taken  to  determine  the  effect  of  marker  density  on 
simulation  convergence.  The  CTH  marker  simulation  results  using  6  markers  per 
cell  were  nearly  identical  to  the  results  of  the  simulations  using  3  markers  per  cell. 
Consequently,  it  was  decided  that  for  the  remainder  of  the  study,  3  markers  per  cell 
would  be  used  for  all  CTH  marker  simulations  and  those  results  are  reported  herein. 

As  previously  described,  ALE  simulations  were  performed  to  compare  to  the  CTH 
Euler  and  marker  simulations.  The  ALE  simulations  performed  in  the  study  used 
the  same  mesh  size  as  the  CTH  simulations,  with  10  elements  across  the  radius  of 
the  impactor.  In  these  simulations,  the  finite  element  mesh  was  conformal  to  the 
outer  shape  of  the  impactor  and  mesh  relaxation  was  enabled  on  the  interior  of  the 
impactor  body  to  reduce  element  distortion  in  regions  of  high  deformation. 

5.  Results  for  OFHC  Copper  Impactor 

Numerical  simulations  were  performed  to  reproduce  the  conditions  of  the  OFHC 
copper  experiments  shown  in  Table  1.  Separate  sets  of  simulations  were  performed 
using  the  3  constitutive  models  described  earlier.  The  final  impactor  shapes  from  the 
simulations  using  the  Johnson-Cook,  Zerilli-Armstrong,  and  Steinberg-Lund  con- 
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stitutive  models  are  shown  in  Figs.  1,  2,  and  3,  respectively.  Similarly,  the  tabulated 
results  for  the  3  sets  of  OFHC  impactor  simulations  are  summarized  in  Tables  2,  3, 
and  4. 

The  final  impactor  shapes  from  the  simulations  using  the  Johnson-Cook  constitutive 
model  are  compared  to  the  experiment  in  Fig.  1  in  which  the  Euler  simulation  is  on 
the  left,  the  marker  simulations  are  in  the  middle,  and  the  ALE  simulation  in  on  the 
right.  The  simulation  results  for  Vs=130  m/s  are  shown  on  the  top  of  the  figure  and 
the  results  for  Vs=190  m/s  are  on  the  bottom.  This  figure  demonstrates  that  all  of 
the  computational  methods  produce  similar  characteristic  results:  1)  the  impact  face 
of  the  body  has  a  larger  diameter  than  the  experiment  with  the  marker  simulations 
producing  a  greater  diameter  than  the  Euler  and  ALE  simulations,  2)  the  diameter 
of  most  of  the  remainder  of  the  body  is  less  than  that  of  the  experiment,  and  3)  the 
final  length  of  the  body  in  the  simulations  is  greater  than  that  of  the  experiment. 

The  results  for  the  OFHC  copper  simulations  using  the  Johnson-Cook  constitutive 
model  are  summarized  in  Table  2.  The  experimental  measurements  of  the  deformed 
copper  penetrators  from  Table  1  are  also  included  in  this  table  for  comparison. 
The  values  of  the  average  error  metric  (A)  for  each  simulation  in  Table  2  show 
the  average  error  for  the  marker  simulation  at  each  striking  velocity  is  greater  than 
those  of  the  Euler  and  ALE  simulations.  Overall,  these  differences  are  the  result 
of  the  impact  face  of  the  marker  simulations  having  a  greater  impact  face  diameter 
than  the  Euler  and  ALE  computational  methods. 

In  conducting  similar  numerical  simulations  of  these  impact  events,  Johnson  and 
Holmquist  made  the  observation  that  simulation  results  for  the  OFHC  copper  im¬ 
pactor  were  improved  by  using  the  Zerilli-Armstrong  constitutive  model  in  place 
of  the  Johnson-Cook  model.9  Accordingly,  it  was  decided  to  repeat  the  previously 
mentioned  simulations  for  the  current  study  using  this  model.  Fig.  2  illustrates  the 
final  shapes  of  the  impactor  from  these  simulations  as  compared  to  the  experiments, 
and  Table  3  contains  a  summary  of  the  results.  These  results  confirm  the  findings 
of  Johnson  and  Holmquist  that  the  use  of  the  Zerilli-Armstrong  constitutive  model 
significantly  improved  the  final  impactor  shape  in  the  simulations  as  compared  to 
the  experiments. 
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Vs=130  m/s 


Vs=190  m/s 


Fig.  1  Final  shape  of  OFHC  copper  impactor  for  experiments  and  simulations  with  Johnson- 
Cook  constitutive  model  (Top:  Vs=130  m/s,  Bottom:  Vs=190  m/s) 


Table  2  Simulation  results  for  OFHC  copper  using  Johnson-Cook  constitutive  model 


Method 

vs 

(m/s) 

L 

(cm) 

D 

(cm) 

W 

(cm) 

A 

Euler 

130 

2.044 

1.098 

0.877 

0.059 

Marker 

130 

2.060 

1.231 

0.880 

0.105 

ALE 

130 

2.064 

1.128 

0.886 

0.069 

Experiment 

130 

1.956 

0.991 

0.899 

Euler 

190 

1.692 

1.431 

0.933 

0.059 

Marker 

190 

1.703 

1.604 

0.921 

0.108 

ALE 

190 

1.730 

1.463 

0.941 

0.072 

Experiment 

190 

1.621 

1.351 

1.007 
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Vs=130  m/s 


Vs=190  m/s 


Fig.  2  Final  shape  of  OFHC  copper  impactor  for  experiments  and  simulations  with  Zerilli- 
Armstrong  constitutive  model  (Top:  Vs=130  m/s,  Bottom:  Vs=190  m/s) 


Table  3  Simulation  results  for  OFHC  copper  using  Zerilli-Armstrong  constitutive  model 


Method 

vs 

(m/s) 

L 

(cm) 

D 

(cm) 

W 

(cm) 

A 

Euler 

130 

1.983 

1.021 

0.896 

0.016 

Marker 

130 

1.984 

1.162 

0.903 

0.064 

ALE 

130 

1.995 

1.074 

0.906 

0.037 

Experiment 

130 

1.956 

0.991 

0.899 

Euler 

190 

1.651 

1.291 

0.973 

0.032 

Marker 

190 

1.650 

1.473 

0.975 

0.047 

ALE 

190 

1.677 

1.349 

0.980 

0.021 

Experiment 

190 

1.621 

1.351 

1.007 

9 


Vs=130  m/s 


Vs=190  m/s 


Fig.  3  Final  shape  of  OFHC  copper  impactor  for  experiments  and  simulations  with  Steinberg- 
Lund  constitutive  model  (Top:  Vs=130  m/s,  Bottom:  Vs=190  m/s) 


Table  4  Simulation  results  for  OFHC  copper  using  Steinberg-Lund  constitutive  model 


Method 

vs 

(m/s) 

L 

(cm) 

D 

(cm) 

W 

(cm) 

A 

Euler 

130 

2.003 

1.008 

0.891 

0.017 

Marker 

130 

2.004 

1.157 

0.897 

0.064 

ALE 

130 

2.013 

1.054 

0.904 

0.033 

Experiment 

130 

1.956 

0.991 

0.899 

Euler 

190 

1.676 

1.249 

0.968 

0.160 

Marker 

190 

1.671 

1.433 

0.965 

0.222 

ALE 

190 

1.703 

1.298 

0.982 

0.038 

Experiment 

190 

1.621 

1.351 

1.007 
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Finally,  the  OFHC  copper  simulation  results  using  the  Steinberg-Lund  constitutive 
model  are  presented  in  Fig.  3  and  Table  4.  These  results  were  nearly  identical  to  the 
results  of  the  simulations  using  the  Zerilli-Armstrong  model,  with  the  Steinberg- 
Lund  model  appearing  to  result  in  slightly  greater  impactor  stiffness  than  the  Zerilli- 
Armstrong  model. 

The  simulations  using  the  Zerilli-Armstrong  and  Steinberg-Lund  constitutive  mod¬ 
els  demonstrated  a  similar  characteristic  as  those  in  the  Johnson-Cook  simulations 
in  that  the  marker  method  produced  a  greater  impact  face  diameter  than  the  Euler 
and  ALE  methods.  Overall,  each  of  the  computational  methods  produced  acceptable 
results  as  compared  to  experiment  providing  a  strong  indication  that  the  constitutive 
models  are  correctly  implemented  into  the  CTH  Lagrangian  marker  method. 

6.  Results  for  4340  Steel  (Rc30)  Impactor 

In  the  next  phase  of  the  study,  numerical  simulations  were  performed  to  match  the 
conditions  of  the  4340  steel  impact  experiments  summarized  in  Table  1.  As  be¬ 
fore,  Euler,  marker,  and  ALE  simulations  were  performed  for  each  impact  velocity. 
Only  the  Johnson-Cook  and  Zerilli-Armstrong  constitutive  models  were  used  in  this 
phase  of  the  study.  The  Steinberg-Lund  constitutive  model  was  not  used  here  be¬ 
cause  parameters  for  4340  steel  with  a  hardness  of  Rc30  were  not  readily  available 
for  this  model.  The  final  shapes  of  the  impactors  from  the  Johnson-Cook  simula¬ 
tions  are  compared  to  the  experimental  shapes  in  Fig.  4,  and  the  tabulated  results 
are  shown  in  Table  5.  Similarly,  the  results  for  the  Zerilli-Armstrong  simulations 
are  provided  in  Fig.  5  and  Table  6. 

The  4340  steel  simulations  produced  results  that  are  consistent  with  the  characteris¬ 
tics  of  the  OFHC  copper  simulations.  The  Zerilli-Armstrong  constitutive  model  pro¬ 
duced  final  impactor  shapes  that  better  matched  the  experiments  than  the  Johnson- 
Cook  model.  Regarding  the  average  error  metric  (A)  of  the  4340  steel  simulations, 
no  single  computational  method  produced  the  minimum  error  across  the  range  of 
impact  velocities  and  constitutive  models  considered.  The  overall  consistency  in  the 
final  impactor  shapes  between  the  3  computational  methods  provides  confidence  in 
the  proper  implementation  and  functioning  of  the  constitutive  models  in  the  marker 
code. 
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V  =208  m/s 


V  =343  m/s 


Fig.  4  Final  shape  of  4340  steel  impactor  for  experiments  and  simulations  with  Johnson-Cook 
constitutive  model  (Top:  Vs=208  m/s,  Bottom:  Vs=343  m/s) 


Table  5  Simulation  results  for  4340  steel  using  Johnson-Cook  constitutive  model 


Method 

vs 

(m/s) 

L 

(cm) 

D 

(cm) 

W 

(cm) 

A 

Euler 

208 

2.253 

1.078 

0.802 

0.013 

Marker 

208 

2.268 

1.199 

0.792 

0.057 

ALE 

208 

2.257 

1.119 

0.809 

0.027 

Experiment 

208 

2.253 

1.062 

0.831 

Euler 

343 

0.604 

1.147 

0.941 

0.025 

Marker 

343 

0.589 

1.288 

0.887 

0.072 

ALE 

343 

0.566 

1.266 

1.006 

0.103 

Experiment 

343 

0.631 

1.147 

0.912 
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V  =208  m/s 


Fig.  5  Final  shape  of  4340  steel  impactor  for  experiments  and  simulations  with  Zerilli- 
Armstrong  constitutive  model  (Top:  Vs=208  m/s,  Bottom:  Vs=343  m/s) 


Table  6  Simulation  results  for  4340  steel  using  Zerilli-Armstrong  constitutive  model 


Method 

vs 

(m/s) 

L 

(cm) 

D 

(cm) 

W 

(cm) 

A 

Euler 

208 

2.263 

0.978 

0.828 

0.029 

Marker 

208 

2.277 

1.109 

0.818 

0.023 

ALE 

208 

2.263 

1.004 

0.840 

0.023 

Experiment 

208 

2.253 

1.062 

0.831 

Euler 

343 

0.628 

1.074 

0.927 

0.028 

Marker 

343 

0.619 

1.228 

0.906 

0.032 

ALE 

343 

0.626 

1.139 

0.975 

0.066 

Experiment 

343 

0.631 

1.147 

0.912 
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For  each  of  the  impactor  materials  and  constitutive  models  considered,  the  marker 
simulations  produced  a  noticeably  greater  impact  face  diameter  than  the  experi¬ 
ments  and  the  other  simulation  methods  for  both  impact  velocities  considered.  This 
increased  impact  face  diameter  is  a  result  of  material  failure  in  the  marker  simula¬ 
tions  that  does  not  occur  as  extensively  in  the  other  computational  methods.  This 
effect  is  illustrated  in  Fig.  6  which  shows  the  distribution  of  damage  in  the  final 
shape  of  the  4340  steel  impactor  for  the  simulations  using  the  Zerilli-Armstrong 
constitutive  model  for  Vs=343  m/s.  The  damage  in  the  Euler  simulation  is  well  be¬ 
low  1 .0 — the  value  at  which  the  material  is  considered  to  be  failed  and  can  no  longer 
support  shear.  The  damage  in  the  ALE  simulation  is  almost  entirely  below  the  fail¬ 
ure  level,  with  only  a  small  portion  of  the  impact  face  exhibiting  failure.  However, 
in  the  marker  simulation,  nearly  the  entire  impact  face  is  failed,  resulting  in  a  larger 
amount  of  radial  deformation  as  compared  to  the  Euler  and  ALE  simulations.  Fur¬ 
thermore,  failure  in  the  marker  simulation  also  occurs  away  from  the  impact  face, 
which  is  not  observed  in  the  Euler  and  ALE  simulations. 

Damage 

1.000e+00 

0.75 

0.5 
^0.25 


Fig.  6  Damage  in  final  impactor  shape  for  4340  steel  simulations  using  Zerilli-Armstrong 
constitutive  model  (Vs=343  m/s) 


Euler  Marker  ALE 


Examination  of  the  early  deformation  behavior  of  the  impactor  provides  insight  into 
the  cause  of  the  increased  extent  of  failure  in  the  marker  simulation  as  compared 
to  the  Euler  and  ALE  simulations.  The  length  of  the  4340  steel  impactor  used  in 
the  Vs=343  m/s  experiment  was  8.13  mm.  For  a  wave  speed  of  4.6  mm///s,  the 
time  for  a  shock  wave  to  propagate  from  the  impact  face  to  the  rear  of  the  body 
is  approximately  1.77  /is.  Fig.  7  provides  images  of  the  deformed  body  at  a  time 
of  2.9  fis  after  impact — sufficient  time  for  a  shock  wave  to  travel  to  the  rear  of  the 
body  and  a  rarefaction  wave  to  begin  traveling  back  toward  the  impact  face. 
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Fig.  7  Marker  simulation  results  for  4340  steel  impactor  at  t=2.9  (is  after  impact  (Vs=343  m/s, 
Zerilli- Armstrong  constitutive  model) 


The  image  on  the  left  side  of  Fig.  7  shows  regions  of  localized  damage  forming  as  a 
result  of  the  shock  wave  propagation  through  the  body.  The  image  on  the  right  of  the 
figure  shows  regions  of  high  equivalent  plastic  strain  occurring  at  the  failure  sites 
after  the  initial  shock  wave  has  propagated  through  the  body.  These  images  demon¬ 
strate  a  dependence  of  the  Johnson-Cook  fracture  model  behavior  on  the  computed 
pressure  field.  To  test  this  hypothesis,  an  additional  marker  simulation  was  per¬ 
formed  using  the  Zerilli-Armstrong  constitutive  model  in  which  the  pressure-stress 
coefficient  in  the  Johnson-Cook  fracture  model  (D3  in  Eq.  11)  was  set  to  zero  to 
eliminate  the  pressure-dependent  portion  of  the  fracture  model.  Fig.  8  is  a  plot  of 
the  final  shape  of  the  impactor  from  this  simulation  which  shows  only  minimal  fail¬ 
ure  at  the  impact  surface,  and  no  failure  interior  to  the  body.  This  result  confirms 
the  assumption  that  the  excessive  failure  in  the  marker  simulations  was  a  result  of 
the  coupling  between  the  pressure  field  and  the  Johnson-Cook  fracture  model. 

7.  Conclusion 

A  computational  study  was  performed  to  evaluate  the  functionality  of  the  Lagrangian 
marker  method  recently  implemented  into  the  CTH  shock  physics  code.  Simula¬ 
tions  were  performed  to  reproduce  the  conditions  of  selected  Taylor  impact  exper¬ 
iments  to  verify  the  proper  implementation  of  3  commonly  used  metal  constitu¬ 
tive  models  in  the  CTH  Lagrangian  marker  method.  Additional  simulations  were 
performed  using  the  legacy  CTH  Eulerian  method  and  an  arbitrary  Lagrangian- 
Eulerian  code  to  compare  material  response  and  isolate  the  influences  of  computa¬ 
tional  method  and  constitutive  model  on  the  overall  results. 
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Fig.  8  Damage  in  final  impactor  shape  for  4340  steel  marker  simulation  using  Zerilli- 
Armstrong  constitutive  model  with  D3=0  (Vs=343  m/s) 


The  study  results  demonstrated  overall  consistent  material  response  between  the 
3  computational  methods  for  each  of  the  constitutive  models  evaluated,  indicating 
correct  implementation  of  the  constitutive  models  in  the  marker  method.  Simula¬ 
tions  using  the  CTH  Lagrangian  marker  method  consistently  resulted  in  greater  im¬ 
pact  face  deformation  than  the  other  2  methods.  Investigation  of  the  data  revealed 
that  the  marker  simulations  had  a  greater  degree  of  material  failure  than  the  Euler 
and  ALE  simulations.  This  material  failure  in  the  marker  simulations  is  caused  by 
the  initial  shock  wave  propagation  through  the  material  and  is  therefore  related  to 
the  relationship  between  the  equation  of  state  and  the  fracture  model.  This  coupling 
between  the  computation  of  the  pressure  field  and  the  fracture  model  requires  fur¬ 
ther  investigation  before  the  CTH  Lagrangian  marker  method  can  be  applied  to  a 
wider  range  of  problems  in  terminal  ballistics. 


16 


8.  References 


1.  McGlaun  J,  Thompson  S.  CTH:  A  three-dimensional  shock  wave  physics 
code.  International  Journal  of  Impact  Engineering.  1990;10:351-360. 

2.  Kimsey  K,  Schraml  S,  Hertel  E.  Scalable  computations  in  penetration  me¬ 
chanics.  Advances  in  Engineering  Software.  1998;29(3-6):209-215. 

3.  Bolling  D.  Evaluation  of  the  material  point  method  within  CTH  to  model  2- 
dimensional  plate  impact  problems.  U.S.  Army  Research  Laboratory;  2014 
September.  Report  No.:  ARL-MR-0880. 

4.  Schumacher  S,  Ruggirello  K.  CTH  reference  manual:  Marker  technologies. 
Albuquerque  (NM):  Sandia  National  Laboratory;  2013  March.  Report  No.: 
SAND2013-1675. 

5.  House  J.  Taylor  impact  testing.  Eglin  AFB  (FL):  Air  Force  Armament  Labora¬ 
tory;  1989  January.  Report  No.:  AFATL-TR-89-41. 

6.  Nichols  A.  Users  manual  for  ALE3D:  An  arbitrary  Lagrange/Eulerian  2D  and 
3D  code  system.  Livermore  (CA):  Lawrence  Livermore  National  Laboratory; 
2010  Jun.  Report  No.:  LLNL-SM-433954. 

7.  Johnson  G,  Cook  W.  A  constitutive  model  and  data  for  metals  subjected 
to  large  strains,  high  strain  rates  and  high  temperatures.  In:  Ballistics;  1983 
Aprl9-21;  The  Hague,  The  Netherlands,  p.  541-547. 

8.  Johnson  G,  Holmquist  T.  Test  data  and  computational  strength  and  fracture 
model  constants  for  23  materials  subjected  to  large  strains,  high  strain  rates, 
and  high  temperatures.  Los  Alamos  (NM):  Los  Alamos  National  Laboratory; 
1989.  Report  No.:  LA-1 1463-MS. 

9.  Johnson  G,  Holmquist  T.  Evaluation  of  cylinder-impact  test  data  for  constitu¬ 
tive  model  constants.  Journal  of  Applied  Physics.  1988;64(8):3901-3910. 

10.  Zerilli  F,  Armstrong  R.  Dislocation-mechanics-based  constitutive  relations  for 
material  dynamics  calculations.  Journal  of  Applied  Physics.  1987 ;6 1(5):  1 8 1 6— 
1825. 

1 1 .  Gray  G,  Chen  S,  Wright  W,  Lopez  M.  Constitutive  equations  for  annealed  met¬ 
als  under  compression  at  high  strain  rates  and  high  temperatures.  Los  Alamos 
(NM):  Los  Alamos  National  Laboratory;  1994.  Report  No.:  LA-12669-MS. 


17 


12.  Steinberg  D,  Lund  C.  A  constitutive  model  for  strain  rates  from  10  4  to 
106  s~4.  Journal  of  Applied  Physics.  1989;65(4):  1528—1533. 

13.  Steinberg  D,  Cochran  S,  Guinan  M.  A  constitutive  model  for  metals  applicable 
at  high-strain  rate.  Journal  of  Applied  Physics.  1980;5 1(3):  1489-1504. 

14.  Steinberg  D.  Equation  of  state  and  strength  properties  of  selected  metals.  Liver¬ 
more  (CA):  Lawrence  Livermore  National  Laboratory;  1991  Lebruary.  Report 
No.:  UCRL-MA- 106439. 


18 


1  DEFENSE  TECHNICAL 
(PDF)  INFORMATION  CTR 

DTIC  OCA 

2  DIRECTOR 

(PDF)  US  ARMY  RESEARCH  LAB 
RDRL  CIO  LL 

IMAL  HRA  MAIL  &  RECORDS  MGMT 

1  GOVT  PRINTG  OFC 
(PDF)  A  MALHOTRA 

3  NAVSEA  DAHLGREN 
(PDF)  W  CHEPREN 

CDYKA 
M  HOPSON 


9  SANDIA  NATIONAL  LABORATORIES 

(PDF)  S  ATTAWAY 

A  BRUNDAGE 
A  GULLERUD 
E  HERTEL 
J  KORBIN 
D  LITTLEWOOD 
KMISH 

S  SCHUMACHER 
S  SILLING 


2  LLNL 
(PDF)  D  FAUX 

R  MCCALLEN 


B  AYDELOTTE 
D  MALLICK 
C  MEYER 
D  SCHEFFLER 
B  SCHUSTER 
RDRL  WMM  B 
BLOVE 
RDRL  WMP  A 
S  BILYK 
RDRL  WMP  B 
K  ZIEGLER 
RDRL  WMP  C 
R  BECKER 
B  LEAVY 
S  SEGLETES 
RDRL  WMP  D 
W  CLARK 
R  DONEY 
D  KLEPONIS 
H  MEYER 
F  MURPHY 
C  RANDOW 
I RUNYEON 
G  VUNNI 
M  ZELLNER 
RDRL  WMP  E 
B  CHAMISH 
D  HORNBAKER 
J  HOUSKAMP 
RDRL  WMP  G 
R  BANTON 
I  STEWART 


4  SOUTHWEST  RESEARCH  INSTITUTE 

(PDF)  C  ANDERSON 
G  JOHNSON 
T  HOLMQUIST 
J  WALKER 


1  UNIV  OF  ALABAMA  BIRMINGHAM 

(PDF)  D  LITTLEFIELD 


2  DSTL 
(PDF)  A  MEARNS 
J  CORDELL 


ABERDEEN  PROVING  GROUND 

26  DIR  USARL 
(PDF)  RDRL  CIH  S 

J  CAZAMIAS 
RDRL  WML  H 


19 


Intentionally  left  blank. 


20 


